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1. Introduction 

The importance of eolian erosion and deposition on Mars is now vm- 
questioned. The studies undertaken during the past four years under this 
project have provided some insight into the nature of eolian processes 
from a viewpoint which has been nearly neglected in the past. We believed 
that in order to interpret the eolian land forms of Mars, it is necessary 
to have a sufficient understanding of the physics of eolian transport and 
deposition on Earth to develop scaling laws which will allow us to scale 
processes and forms to the vastly different surface and atmospheric regime 
of Mars. In peurticular, would the eolian landforms found on Earth be 
present under the thinner atmosphere and stronger wind regime, and if 
present, how would the size and shape of these land forms compare with 
those on Earth? 

During the first three years of this research one of us (A.H.) 
conducted field and aerial photography surveys of some eolian-dominated 
landforms on earth (in the United States and Peru). During the past year 
the emphasis has been on wind tunnel studies of the flow field around 
models of these eolian forms to attempt to reproduce a more extensive map 
of the flow. This effort has given us information useful in elucidating 
the mechanics determining the form and scale of eolian bedforms and the 
Interactions between topography and regional winds which result in sand 
deposits. The shift in research emphasis reflects the need for quantitative 
verification and extension of hypotheses generated during the field studies. 
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2. Areas of Activities 

The wind tunnel studies have been directed primarily in the following 
five areas: 

I. Simulation of the atmospheric boundary layer in a wind tunnel, 

II —Velocity profile measurements around different models in the desert 
boundary layer; estimating shear stress distributions on the model 
surfaces, 

III —Developing appropriate flow visulization techniques for better 
understanding of the flow around the different models, 

IV- Streamline mapping using tuft photographs, 

V— Roughness contrast experiments. 
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3. Simulation of Atmospheric Boundary Layer 
in a Wind Tunnel 

Techniques have been developed for simulating the lower 100 meters 
or so of a neutral atmosphere. The initial effort was directed at adapting 
these techniques to flow over a desert. This Involved simulating an 
appropriate velocity profile and ground roiighness. 

In the wind tunnel used, the incoming air is initially forced through 
a contraction (C.P. figure l), producing a more uniform velocity profile 
with turbulence levels at approximately 0.1J5. It then passes through two 
grids (#1 and H2) designed to produce high turbulence levels similar in 
magnitude to those observed in the atmosphere. A final grid of variable- 
spaced bars (#3) produces the desired logarithmic velocity profile, correspond- 
ing to the desired z© [l]. The models are embedded into the roughened 
floor . 

For a neutrally stable atmosphere, the mean velocity profile near 
the ground (the first 100 meters or so, where Coriolis effects are negligible) 
can be written as: 

U (z) = — • n ~ 

K Zo 

where u- is the friction velocity (= '^o, x is the shear stress at the 

* / 0 

'f V 

surface euid p is the air density), k is the von Karman constant (; O.l*), 
and z^ is the roughness height [2]. 

A convenient scale factor for our facilities has been found to be 
315 :!• Thus, different dunes have been scaled down by that factor; an 
atmospheric boundary layer of height 190 m. is roughly represented by our 
tunnel height, 60 cm., and a natural roughness of 7.2 cm. is simulated in 
our tunnel by 0,023 cm. Figure 2 is a semi-logarithmic plot of the mean 
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velocity profiles at different distances (x) downstream from the logarithmic 
rod grid (rod diameter, d, is I .9 cm.}. The profiles follow the given 
logarithmic velocity profile formula. Figure 3 shows the turbulence level 
profiles. This simulation of the atmospheric boundary layer is quite 
satisfactory and ccxnpares well with naturally grown boundary layers in wind 
tunnels [l, 3, 5]. 
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4.. Velocity Profile Meaauremanta Around 
Different Models In the Desert Boundary Layer 

Velocity profile measurements were made around four different topographic 
forms In the simuilated desert flow field of the wind tunnel. The first 
three models are scale models of harchan dunes that exist in an area near 
the Salton Sea, California; the fourth model is a conical "hill” with a 
slope of 33 degrees. Detailed mean velocity profiles at different 
locations on the cone and on cne of the harchan dvmes were measured. 

Figure U shows the locations at which profiles were measured on the dune, 
called harchan #1 for convenience, a contotir map of harchan #1 is shown in 
figure 5. Profiles were taken at the intersections of lines A, A’ (the 
dime's crest), B, and C with lines 1 through 13; additional profiles were 
measured at various locations in front of and behind the dune. A less 
extensive set profiles was measured around the two other harchan dune 
models. 

The profiles were measured at a constant reference tunnel speed 

10 m/s; the effect of changing the reference speed, U , was studied, and 

n 

it will he shown later that speed changes in the range 5-10 m/s have no 
effect on the non-dimensionalized profiles. Only the component of the 
mean velocity profile parallel to the flow direction in the wind tunnel 
was measured, although lateral velocities could he estimated from the long- 
itudinal component and the streamline maps resulting from flow visualization 
measurements. 

The models were inserted at a distance 287 cm downstream from the 
rod grid. Reference txinnel speed was monitored hy a Mariam manometer 
model 3^*FBZ connected to two static pressure holes in the tunnel walls. 

Mean velocity profiles and turbulence intensities were measured using 
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two types of hotwire probes, a Disa type 55F11* 5 micron, platinum- plated, 
tungsten, stralghtwire, with a sensitive length of 1.23 mm, and a Disa type 
55P12 probe, which facilitates measurements close to the windward 
surfaces of the topographic forms. A Disa type 33D01 constant tenperature 
anemometer was used in conjunction with a Disa type 3^D01 linearizer and 
a Disa type 55D29 auxiliary unit. An overheat ratio of 0.6 and a 
linearizer exponent setting of 2.2 were normeLlly used, producing an output 
linearly proportional to the velocity fluctuations within a set 
frequency range from D.C»to 5 KHz. 

Figure 6 shows a composite dieigram of selected velocity profiles 
taken around barchan #1. Figure 7 shows a composite of velocity profiles 
around the cone. In each diagram, an undisturbed velocity profile, U, 
measured 10 cm upwind of the models, is included for comparison. The 
The comi>OBite profiles include only a few of the profiles measured and 
were selected because of the significant flow trends they show. The 
profi'’es not shown in the composite diagrams are being subjected to 
different types of analysis, to be described later. 

From a stud; of figure 6, it is possible to deduce a few flow 
characteristics along line 7. (Refer to figure 1^). Velocity increases 
as the flow travels up the windward face of the dune toward the crest. 

All of the profiles rejoin the undisturbed logarithmic velocity profiles 
at a height of approximately 12 cm above the ground; thus the dune 
significantly affects the flow up to a height 6 times its own height. A 
"kink" in all of the curves occurs at a height 10 cm above the ground. The 
kink was a repeatable phenomenon and unexpected; we are unable to explain 
its appearance or significance at the present time. In the lee of the 
dune, measurements are difficult due to the inability of a hot wire to 
follow reverse flows [l7l. However, the reverse flow region extends 
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only to about 2 cm above the ground, and from the profiles one could 
estimate the location of the "dead zone" behind the dune. Similar 
observations can be made about the flow around the cone [6, 7] (figure T). 
Additional composite diagrams have been carefully studied for flow trends. 

Additional experiments were performed in studying the flow around 
Barchan #1. As mentioned previously, the tunnel reference speed was 

I changed to investigate any possible effects on the velocity profiles. 

f 

I Profiles at several points were measured for tunnel speeds. Up, of 5» 

6, 7, 8 and 10 m/s. It was found that if the velocity profiles were 
nondimens ionali zed using the appropriate reference speed, the profiles 
would collapse with a percentage: error of about UJK, on the average. 

Figure 8 shows sample velocity profiles taken (at point 7A) at reference 
speeds of 5 and 10 m/s; the two profiles were nondimens ionali zed by 
dividing the velocities by the appropriate reference speed. It can be 
seen that the profiles closely coincide. 

Another experiment was performed to determine the necessity of using 
the simulated boundary layer for measuring barchan profiles. The grids 
and roughness element were removed from the tunnel , and profiles were 
taken at several points on the dune. Figure 9 shows a typical comparison 
between the free-stream profile at point 7A and the logarithmic profile 
at the same location. It is clear from the graph that the curves are 
entirely dissimilar. Therefore, it was established that the simulated 
environment was necessary if meaningful velocity profiles around the 
barchan dune were to be produced. 

A third experiment was conducted to investigate the effects of change 
of wind direction on the profiles. The dune was rotated through several 
known angles, producing the effect of a change in wind direction, and, at 
a reference tunnel speed of 10 m/s, velocity profiles were taken at several 
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points for each angle. Analysis of the results of this experiment is still 
not completed, alttough a sheltering effect has been noted at certain 
points when successive angles of rotation move the point more in the lee 
of the dune's bulk. Figures 10 and 11 show the sheltering effect 
describedi at position 13A' (on the crest of the dune), as the dune is 
rotated from 30° to -30° the velocities close to the surface decrease as 
the point 13A' is moved to the lee of the main part of the dune, 
sheltering the point from the wind. At 2A'(2 crest) a similar effect is 
evident as the dune is rotated; in this case the velocities decrease 
as the dune is rotated from -30° to 20° because 2A'is on the opposite 
wing of the dune. As z increases, the angular effect ceases, as the 
profiles join the undisturbed logarithmic profile. 
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5. Developing Appropriate Flow Visualization 
Techniques for Better Understanding of 
the Flow Around the Different Models 

Hotwire measurements made around the models provide accurate Information 
about the relative speeds of the flow at different points In space. The 
ability of the hotwire to indicate direction of flow is limited and 
involves a tedious process. The hotwire has an angular cosine sensitivity 
to the flow 6ind It may be rotated to find the direction of flow. However, 
it must be maneuvered back to the position of interest after it is rotated 
through the desired angle, and the prongs on which the wire is mounted can 
cause some Interference with the actual flow if they are blocking the 
sensitive part of the wire. 

An alternate method of determinifit flow direction is flow visual- 
ization. There are many methods of flow visualization used for deter- 
mining streamlines around models; three principal methods were investigated 
for this project. The first, and moat suc'it’ssful, method used was tufting. 
Thin strips of tissue paper were mounted at various points on and around 
the dune, carefully glued in place to have no directional bias initially. 
Under wind action, the tufts aligned themselves in the direction of the 
wind. Photographs were taken at different exposures to determine 
"instantaneous” wind direction as well eus "average" wind direction, and 
the latter were used in the construction oi streamline maps. A sample 
streamline map is included in the next section. Figure 12 (a) shows a 
sample tuft photograph taken at 1/2 second exposure; 12 (b) was taken at 
an exposure time of 2 seconds. 

The second method of flow visualization employed was smoke injection. 

A smoke generator, the schematic design of which is reproduced in figure 13, 
produced a Jet of smoke, initially laminar, which was illuminated on one 
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plane by a laser light. The resulting flow could be seen in that plane 
clearly with the eye and sonevhat less clearly in the photographs. A 
Tri'X film was used and developed with UQF high-speed developer to yield 
an effective ASA number of 3200. The chief drawback of the smoke flow 
visualization was that* although the flow ccUd be seen clearly when the 
only velocity was that of the smoke emerging from the orifice* once the 
tunnel wind was turned on* the smoke dispersed so rapidly that it br'came 
too thin to scatter enough light for an effective photograph. Other 
drawbacks were poor photograph quality and the limitation imposed on the 
planes of illuminated flow due to the inaccessibility of lower planes. 

The third method of flow visualization used was the oil drop method. 
Colored drops of oil were quickly placed on the dune's surface through a 
hypodermic syringe and the wind was turned on. The resulting tracks left 
by the colored oil as it travelled over the dune's surface were to 
have been the streamlines. The oil drop method was successful when a high 
speed Jet of air was aimed at the model. For speeds attainable in the 
wind tunnel (5 - 10 m/s), however, gravitational effects outweighed wind 
effects on the oil drops. 

Thus* tuft flow visualization was ultimately the most successful 
streamline indicator employed of the three flow visualization methods. 
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6. Streamline Mapping Using Tuft Photographs Estimate 
for Shear Strcsa PietributlonB on the Surface 

Estlfflates of erosion and deposition rates of sand on the surTarp of 
a barchan dune are of prinary importance in any analysis to explain the 
self-preserving nature of the dune's shape. Both streamline mapping and 
detailed mapping of the wind velocities close to the dune's surface are 
ci*ucial factors in determining erosion and deposition rates. Using the 
tufting method of flow visualization mentioned in section 5 and the near 
surface velocity profiles described in section U, It is possible to 
construct streamline maps for the barchan dune models. 

Figure lU a map of the near surface streamlines for barchan #1 when 
the tunnel reference speed is 10 m/s. Streamlines on the upwind surface 
of the dune are the result ^’f analysis of tuft photographs. Several 
photographs are used to construct one map; average tuft direction is deter- 
mined for each tuft from the photographs, and the resulting direction 
vectors are plotted on a diagram of the dune. The vectors are extrapolated 
end Joined to form streamlines. Allen [8] used flow visualization methods 
in a water flume to find streamlines for a dune similar to the barchan 
represented in figure lU. Allen's streamlines (dashed lines) for the wind 
fli'' in the lee of the dune are ad led to figure lit to give a complete 
picture of the flow around the dune. 

Quantitative data of the near surface velocity profiles obtained 
(section U) by the hotwire probe provides a good indication of direction 
of wind flow over a barchan dune. Streamline mapping using near velocities 
is still in progress. Measurements as close as 2.5 mm to the dune's surface 
have been made at approximately 50 different points on the model. Sand 
transport occurs chiefly through the mechanism of saltation (9, 10, 11, 
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12, 13, 1^], a process which ejects the grains into a trajectory to an 
average height above the dune's surface. Thus, the near surface velocities 
on the dune are influential in determining the path of the sand grains. 

Detailed calculations relating the near surface velocities and slope angles 
of the dune to the resulting grain paths and, thus, to the total movement 
of sand around the barchan dune, are in progress. Further experimentation 
with tufting to determine streamlines will be carried out to aid in the 
quantitative analysis of the flow. 

The use of shear stress probes, Disa-type 55H^7, has been investigated 
for determining the shear stress at the surface of the dune. Shear stress 
measurements will be compared to near surface velocity measurements to 
determine a possible correlation between them. The principle source of 
error in the use of the shear stress probe has been found to be in its 
calibration process. A sample calibration curve is given in figure 15. 

The probe is calibrated to give a value of shear stress vs. voltage output 
in the following manner: the probe is mounted on a flat surface, perpendicular 

to the oncoming wind; a pitot tube is located close to the probe and is 
moved vertically. For each of several tunnel speeds, a velocity profile 
is measured with the pitot tube. The output voltage from the shear stress 
probe is lineeu’ized and is constant for a constant reference tunnel speed. 

From the velocity profiles, values of the friction velocity, u^,, and hence 
the shear stress x at the surface, are obtained for each reference tunnel 

d 

speed. Values of are plotted against the corresponding values of 
voltage obtained from the shear stress probe. The scatter of the points 
on the calibration curve is a result of the uncertainty inherent xn 
determining from the velocity profiles. 
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7. Roughness Contrast Experiments 

When wind conditions responsible for preserving isolated barchan dunes 
change, a different sand form results. According to Bagnold [9], strong 
winds emerging from a direction other than that of the prevailing gentle 
wind may be the cause of a transition from barchan dune formations to 
longitudinal, or seif, dunes, whio.i form long chains parallel co the wind 
direction. A series of wind tunnel experiments were conducted to determine 
the flow pattern around longitudinal sand forms, which were represented in 
the tunnel by 6 m long strips of roxighness elements with different Zq 
values. Similar experiments had been conducted in the field of agriculture 
to determine the effects of contrast of roughness on soil erosion by 
wind; soil erosion resulting when the contrast occurred perpendicular 
to the wind direction was compared to the erosion when the contrast was 
parallel to the prevailing wind [15, I6]. 

Several combinations of rotighness elements were studied in the wind 
tunnel, using flow visualization technit^ues and pitot tube measurements. 
Stones with an average diameter of l.U cm were aligned in 20 cm wide strips 
on either side of a "smooth” element. The smooth strip varied from cardboard 
to small pebbles (average d = O.U5 cm). Velocity profiles were measured 
at several locations along a line perpendicular to the wind direction. 

An alternate experiment consisted of replacing the smooth middle layer 
with the l.k cm diameter stones, "sandwiched" between two smooth cardboard 
strips. Profiles were measured with the stones aligned both parallel and 
at a slight angle to che wind direction. Figure 16 illustrates two contrast 
arrangements which were of particular interest. Figure I6 (a) shows the 
smooth- rough-smooth combination, and l6 (b) shows the alternate rough-smooth- 
rough arrangement. In both cases, the "smooth" element was cardboard; the 
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”rough" element consisted of the l.U cm diameter stones. Velocity profiles 
were measured at 5 different locations on a line perpendicular to the wind 
direction. 

Figures 17 and l8 are graphs of the profiles measured in case (a) 
above. A trend is evident in the graphs; at z £ 8 cm, the flow is at its 
highest speed at position 5 and progressively decreases in speed as it 
approaches position 1, at the center of the atones. This trend indicates 
a crossflow of wind from positions 1 to 5, which is strongest near the 
floor and weakens with increasing z. At z >_ 8 cm, the profiles converge 
and then reverse slightly (at z = i6 cm), indicating a weak crossflow from 
positions 5 to 1. 

Figures 19 and 20 are graphs of the profiles measured in case (b), 
which also show a crossflow at z 8 cm away from the stones and toward 
the smooth area. The same weak crossflow as for case (a) occurs towards 
the stones at z = IT cm. In cases (a) and (b), the semilog graphs clearly 

show the change in effective Lq as the position change. 

Tuft photographs taken correlate well with the profile measurements, 
indicating a deflection away from the stones (and, thus, a crossflow away 
from them) at low z values , and a slight deflection toward the stones at 
higher elevations. 
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8, Concluding Renarks 

Thw- desired desert toxindary layer, vith a Zq value of 0.23 nun, was 
sLuulated in the wind tunnel and compares satisfactorally with naturally 
grown boundary layers. 

Velocity profiles around four models were measured and significant 
flow trends were observed. Experiments were conducted to determine the 
effects of change of wind speed and direction on the flow. 

Three flow visualization methods were investigated; the "tufting" 
method proved to be the most effective in indicating flow trends. 

Streamline mapping of the barchan dune is in progress, using flow 
visualization methods and the results of near velocity profiles. 
Calculations of sand flow around the barchan are in progress , using the 
results of streamline mapping and estimates for the shear stress. 

Experiments performed on elements of contrasting roughness show 
consistent crossflow patterns. 
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Figure Captions 

1: Schematic of the wind tunnel test section 

2: Undisturbed velocity profiles at various downstream locations 

3: Turbulence level at various downstream locations 

4: Sch«natic grid of profile locations (barchan #l) 

5: Contour map of barchan #1 

6; Velocity profiles along line 7 for barchan #1 
7: Velocity profiles around cone 

8: Scaled velocity profiles at 5 and 10 m/s 

9: Comparison between free-stream and logarithmic velocity profiles 

10: Velocity profiles at 13 crest (13A^) with change of wind angle 

11; Velocity profiles at 2 crest (2A^) with change of wind angle 
12: Tuft flow visualization photographs of barchan #1 

13: Schematic of smoke generator 

lU: Streamline map of barchan #1 

15: Shear stress probe calibration curve 

l6: Roughness contrast experiment 

17: Velocity profiles, inner roughness aisle 

l8: Velocity profiles, inner roughness aisle (semi-log plot) 

19: Velocity profiles, outer roughness aisles 

20; Velocity profiles, outer roughness aisles (semi-log plot) 
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Figure H Schematic of the Wind Tunnel Test Section 





Figure 21 Undisturbed Velocity Profiles at Various Oownstream Locations 
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Figure 3! Turbulence Level at Various Downstream Locations 




Figure 4: Scherwtlc Grid of Profile Locations (Barchan #1) 




Figure 5; Contour Mep of Barchan #1 





Figure oi Velocity Profiles along Line 7 for Barchan II 






Figure 71 Velocity Profiles around Cone 







Figure &s Comparison between Free-Stream and Logarithmic Velocity Profiles 






Figure 13. Schematic of Smoke Generator 





Figure ilf; Streamline Map of Barchan #1 
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Figure 1®: Roughness Contrast Experirrient 
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